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Abstract 

Reliable calculations of financial risk require that the fat-tailed nature of prices changes is in- 
cluded in risk measures. To this end, a non-Gaussian approach to financial risk management is 
presented, modeling the power-law tails of the returns distribution in terms of a Student-t (or Tsal- 
lis) distribution. Non-Gaussian closed-form solutions for Value-at-Risk and Expected Shortfall are 
obtained and standard formulae known in the literature under the normality assumption are re- 
covered as a special case. The implications of the approach for risk management are demonstrated 
through an empirical analysis of financial time series from the Italian stock market. Detailed com- 
parison with the results of the widely used procedures of quantitative finance, such as parametric 
normal approach, RiskMetrics methodology and historical simulation, as well as with previous 
findings in the literature, are shown and commented. Particular attention is paid to quantify the 
size of the errors affecting the risk measures obtained according to different methodologies, by 
employing a bootstrap technique. 
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A topic of increasing importance in quantitative finance is the development of reliable 
methods of measuring and controlling financial risks. Among the different sources of risk, 
market risk, which concerns the hazard of losing money due to the fluctuations of the prices 
of those instruments entering a financial portfolio, is particularly relevant. 

In the financial industry today, the most widely used measure to manage market risk is 
Value-at-Risk (VaR) . In short, VaR refers to the maximum potential loss over a given 
period at a certain confidence level and can be used to measure the risk of individual assets 
and portfolios of assets as well. Because of its conceptual simplicity, VaR has become a 
standard component in the methodology of academics and financial practitioners. However, 
as discussed in the literature , VaR suffers from some inconsistencies: first, it can violate 
the sub-additivity rule for portfolio risk, which is a required property for any consistent 
measure of risk, and, secondly, it doesn't quantify the typical loss incurred when the risk 
threshold is exceeded. To overcome the drawbacks of VaR, the Expected Shortfall (or 
Conditional VaR) is introduced, and sometimes used in financial risk management, as a 
more coherent measure of risk. 

Three main approaches are known in the literature and used in practice for calculating 
VaR and Expected Shortfall. The first method, called parametric or analytical, consists in 
assuming some probability distribution function for price changes and calculating the risk 
measures as closed-form solutions. Actually, it is well known that empirical price returns, 
especially in the limit of high frequency, do not follow the Gaussian paradigm and are char- 
acterized by heavier tails and a higher peak than a normal distribution. In order to capture 
the leptokurtic (fat-tailed) nature of price returns, the historical simulation method is often 
used. It employs recent historical data and risk measures are derived from the percentiles of 
the distribution of real data. A third approach consists in Monte Carlo simulations of the 
stochastic dynamics of a given model for stock price returns and in calculating risk measures 
according to Monte Carlo statistics. 

Actually, reliable and possibly fast methods to calculate financial risk are strongly de- 
manded. Inspired by this motivation, the aim of this paper is to present a non-Gaussian 
approach to market risk management and to describe its potentials, as well as limitations, 
in comparison with standard procedures used in financial analysis. To capture the excess 
of kurtosis of empirical data with respect to the normal distribution, the statistics of price 
changes is modeled in terms of a Student-t distribution (also known as Tsallis distribu- 
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tion |5j), which is known to approximate with good accuracy the distribution derived from 
market data at a given time horizon We present, in the spirit of a parametric approach, 

closed-form expressions for the risk measures (VaR and ES) and critically investigate the 
implications of our non-Gaussian analytical solutions on the basis of an empirical analysis of 
financial data. Moreover, we perform detailed comparisons with the results of widely used 
procedures in finance. Particular attention is paid to quantify the size of the errors affecting 
the various risk measures, by employing a bootstrap technique. 



NON-GAUSSIAN RISK MEASURES 



Value-at-Risk, usually denoted as A*, is defined as the maximum potential loss over a 
fixed time horizon At for a given significance level V* (typically 1% or 5%). In terms of 
price changes AS, or, equivalenty, of returns R = AS/ S, VaR can be computed as follows 

i— A* r- A*/S 

V*= dAS P At (AS) = S dR P At (R), (1) 

J —oo J — oo 

where Pa* (AS) and P At (R) are the probability density functions (pdfs) for price changes 
and for returns over a time horizon At, respectively. Actually, VaR represents the stan- 
dard measure used to quantify market risk because it aggregates several risk component 
into a single number. In spite of its conceptual simplicity, VaR shows some drawbacks, as 
mentioned above. 

A quantity that does not suffer of these disadvantages is the so called Expected Shortfall 
(ES) or Conditional VaR (CVaR), E\ defined as 

E* = ^ [ A dAS (-AS) P At (AS) = A f A /S d R (-R) P M (R), (2) 

' J —oo i J —oo 

with V* and A* as in Eq. (0). 

Assuming returns as normally distributed, i.e. R ~ J\f(m,a 2 ), VaR and ES analytical 
expressions reduce to the following closed-form formulae 

A* = - m S + aS V2 erfc' 1 ^*) (3) 

and 

E* = -mSo + ^-i=exp{-[erfc- 1 (2^)] 2 }, (4) 
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FIG. 1: Convergence of VaR (left) and ES (right) Student-t formulae toward Gaussian. 
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where erfc -1 is the inverse of the complementary error function [8|. Note that expressions (jUJ) 
and (jlj) are linear with respect to the spot price Sn. 

However, it is well known in the literature [3, UOl that the normality hypothesis is often 
inadequate for daily returns due to the leptokurtic nature of empirical data. For this reason, 
a better agreement with data is obtained using fat-tailed distributions, such as truncated 
Levy distributions fl Q „ r Student, one, In o r de r to oh— the excess of k urtos,s, 
we model the returns using a Student-t distribution defined as 

1 a u 



I "+ 1 ' 



(5) 



£0/2,1/2) [ a * + {R 

where v £ (1, +oo) is the tail index and 5(^/2,1/2) is the beta function. It is easy to 
verify that, for v > 2, the variance is given by a 2 = a 2 /(u — 2), while, for v > 4, the excess 
kurtosis reduces to k = 6/(z/ — 4). Under this assumption, we obtain closed- form generalized 
expression for VaR and ES given by 



A" 



and 



E* = -mS Q + 



mS + aSo^/v — 2 

aS \fu 



1-A* 



A* 



i — 

2 



(6) 



(7) 



V*B(v/2, 1/2) v - 1 

where A* = 1 , 2 , (2V*) and x , 2 , is the inverse of the incomplete beta function 
As shown in Fig. Q we have checked numerically the convergence of formulae (jHI) and (J7|) 
to the Gaussian results and (jlj) using different values of tail index v (2.75,3.5,4.5,100). 
As expected, the points corresponding to v = 100 are almost coincident with the Gaussian 
predictions, demonstrating that our results correctly recover the Gaussian formulae as a 
special case. 



We observe that each line, corresponding to a fixed v, crosses over the Gaussian one 
for a certain V* . In the light of this observation, we report in Table H] the values of z/ cross 
corresponding to a given V* for both VaR and ES. As can be observed, the growth of z/ cross 
with V* is very rapid for VaR, while for ES and for usually adopted significance values, z/ cross 
keeps in the interval [2.09,2.51]. From this point of view, VaR and ES are quite different 
measures of risk, since the crossover values for the latter are much more stable than those 
associated to the first one. This result can be interpreted as a consequence of ES as a more 
coherent risk measure than VaR. 

TABLE I: Values of v crossover for VaR and ES for different V*. 

V* 1% 2% 3% 4% 5% 

i/ cross (VaR) 2.44 3.21 5.28 32.38 > 100 
i/ cross (ES) 2.09 2.18 2.28 2.38 2.51 



EMPIRICAL ANALYSIS 

The data sets used in our analysis consist of two financial time series, composed of 
N = 1000 daily returns, from the Italian stock market: one is a collection of data from the 
Italian asset Autostrade SpA (from May 15 th 2001 to May 5 th 2005), while the other one 
corresponds to the financial index Mib30 (from March 27 th 2002 to March 13 th 2006). The 
data have been freely downloaded from Yahoo Finance ^l] . Other examples of analysis of 
Italian stock market data can be found in [Rj ]. 

In order to balance between accuracy and computational time, we estimate mean m 
and variance a from empirical moments. With this standard procedure, we measure daily 
volatilities of the order of 1% for both time series: a Autostrade = 1-38% and aMmo = 1-16%. 
Moreover, we find quite negligible means: m Autostrade = 0.12% and triMmo = 0.02%. Using 
the above m and a values, we derive a standardized vector (with zero mean and unit variance) 
r = (r t , . . . , r t -N+i), where Tt-i = (Rt-i — m)/a for i = 0, . . . , N — 1. In order to find the 
best value for the tail parameter u, we look for the argument that minimizes the negative 
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log-likelihood, according to the formula 



v = argmm 



u>2 



N-l 



V-2 



Tt- 



i=0 



(8) 



where the constraint v > 2 prevents the variance to be divergent and i—^ is as in Eq. © 



with m = and a = \Jv — 2. We remark that the beta function B(u/2, 1/2) only admits an 
integral representation and therefore we implemented a numerical algorithm to search for 
the minimum and solve the optimization problem. We measure tail parameter of 2.91 for 
Autostrade SpA and 3.22 for Mib30. These values confirm the strong leptokurtic nature of 
the returns distributions, both for single asset and market index. 

For completeness, it is worth mentioning that other approaches are discussed in the 



literature to model with accuracy the tail exponent o 



Extreme Value Theory 13] and Hill's estimator 



14, 



id- 



the returns cdfs and are based on 



COMPARISON OF DIFFERENT RISK METHODOLOGIES 

In this Section we present a comparison of the results obtained estimating the market 
risk through VaR and ES according to different methodologies. The standard approach is 
based on the normality assumption for the distribution of the returns. As discussed above, 
we limit our analysis to 1000 daily data and we estimate the volatility using the empirical 
second moment (the effect of the mean is negligible). In order to avoid the problem of a 
uniform weight for the returns, RiskMetrics introduces the use of an exponential weighted 
moving average of squared returns according to the formula (?| 

^m| t = E - ™) 2 ' (9) 

i=0 

where A G (0, 1] is a decay factor. The choice of A depends on the time horizon and, for 
At — 1 day, A = 0.94 is the usually adopted value Q. <Tt+i\t represents volatility estimate at 
time t conditional on the realized R. In order to relax standard assumption about the return 
pdf without loosing the advantages coming from a closed-form expression, we presented 
above generalized formulae for VaR and ES based on a Student-t modeling of price returns. 
As a benchmark of all our results, we also quote VaR and ES estimates following a historical 
approach, which is a procedure widely used in the practice. According to this approach, 
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FIG. 2: Bootstrap histograms for tail index v (left) and for RiskMetrics 
volatility proxy er t+ i|t (right) for Autostrade SpA (M = 10 3 bootstrap copies) 
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after ordering the N data in increasing order, we consider the [iVP*] th return R^-p*}) as an 
estimate for VaR and the empirical mean over first [NV*] returns as an estimate for ES. 

At a variance with respect to previous investigations jy, llfij , we also provide 68% confi- 
dence level (CL) intervals associated to the parameters, to estimate VaR and ES dispersion. 
To this extent, we implement a bootstrap technique |l7[. Given the TV measured returns, we 
generate M = 1000 synthetic copies of R, {RJ}, with j = 1, . . . , M, by random sampling 
with replacement according to the probability p = (1/N, . . . , 1/N). For each R* we estimate 
the quantities of interest 9 and we obtain bootstrap central values as follow 

1 - 

d * b =uJ2 6 *r ( 10 ) 



j'=i 



We define the 1 - 2a CL interval as 0J_ a ], with 9* a such that P{6* < 6* a ) = a and 
a = a, 1 — a. 68% CL implies a = 16%. In Fig. 01 Tables HH and lllll we quote results 



according to 9 



-)\_ a —01). In this way, we use the bootstrap approach in order 



to estimate the dispersion of the quantity of interest around the measured value 9. In Fig. |21 
we present, as example, the histogram of bootstrap values of tail index v for Autostrade 
SpA. We observe that the empirical value v = 2.91 is close to the bootstrap central value 



2.96. We also include 68% CL interval (f* 6% = 2.75 and z/g 4% = 3.16), in order to 
quantify the dispersion around v£. 

As a rule of thumb, we consider the bootstrap approach accurate when, given a generic 
parameter, the difference between its empirical value and the bootstrap central value esti- 
mate is close to zero and 68% CL interval is almost symmetric. In our numerical investi- 
gation, we found a systematic non zero bias for 6 = (Tt+i\t- in Fig. EJ crt+i\t — 1-83% while 



TABLE II: Parameters values and bootstrap estimates for the 68% CL intervals. 



m a a t+l \ t v #( 10 ) 

Autostrade 0.12^;°f% 1.38±g;?g% 1.83±g;||% 2.91^° -3.51±g;?|% 
Mib30 0.02±g;g|% 1.16±g;gg% 0.72+°;|% 3.22ig;^ -3.33^;|% 



FIG. 3: VaR A* (upper panel) and ES E* (lower panel) central val- 

ues with 68% CL intervals for Autostrade SpA (left) and for Mib30 (right). 
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1.32%, so we measured a bias of order 0.005. It is worth noticing the positive 



skewness of the histogram. From Table HP it is quite evident the asymmetry of Rn^-p*]). 
Therefore, we can consider the corresponding CL intervals as a first approximation of the 
right ones 

In Fig. El we show VaR and ES central values and 68% CL bars for Autostrade SpA and 
Mib30, corresponding to 1% and 5% significance level and according to the four method- 
ologies previously described. In Tables IIHI we detail all the numerical results. As already 
noted in Ref. |l6j], at 5% significance level Student-t and Normal approaches are substan- 
tially equivalent, but here such a statement sounds more statistically robust, thanks to the 
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TABLE III: Estimated VaR and ES values (mean and 68% CL inter- 
val) for 1% and 5% significance levels from Autostrade SpA and Mib30. 
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bootstrap 68% confidence levels and to the comparison with the historical simulation. We 
note also, from Fig. |3] and Table IHIl that A* and E* central values calculated according to 
RiskMetrics methodology are quite fluctuating and characterized by the largest CL bars. 
The decreasing of V* traduces in a major differentiation of the different approaches. In 
general, we obtain the best agreement between the Student-t approach and the historical 
simulation, both for A* and E*, whereas, as before, the RiskMetrics methodology overes- 
timates or underestimates the results of the historical evaluation and is affected by rather 
large uncertainties. 

CONCLUSIONS 

In this paper we have presented a careful analysis of financial market risk measures in 
terms of a non-Gaussian (Student-like) model for price fluctuations. With the exception of 
Gaussian ones, the derived closed-form parametric formulae are able to capture accurately 
the fat-tailed nature of financial data and are in good agreement with a full historical 
evaluation. We also proposed a bootstrap-based technique to estimate the size of the errors 
affecting the risk measures derived through the different procedures, in order to give a sound 
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statistical meaning to our comparative analysis. 

Possible perspectives concern the extension of our analysis to other time series, different 
financial instruments and underlying distributions. 
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